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In the paper, we present a new kind of function photonic crystals, which refractive index is a 
function of space position. Unlike conventional PCs, which structure grow from two materials, A 
and B, with different dielectric constants ea and eb- Based on Fermat principle, we give the motion 
equations of light in one-dimensional, two-dimensional and three-dimensional function photonic 
crystals. For one-dimensional function photonic crystals, we investigate the dispersion relation, 
band gap structure and transmissivity, and compare them with conventional photonic crystals. By 
O , choosing various refractive index distribution function n(z), we can obtain more wider or more 

■ narrower band gap structure than conventional photonic crystals. 
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1. Introduction 
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Photonic crystals (PCs), proposed by Yablonovitch and John, represent a novel class of optical materials 
which allow to control the flow of electromagnetic radiation or to modify light- matter interaction [1, 2]. 
These artificial structures are characterized by one, two or three-dimensional arrangements of dielectric 
material which lead to the formation of an energy band structure for electromagnetic waves propagating in 



them. One of the most attractive features of photonic crystals is associated with the fact that PCs may 
Oh exhibit frequency ranges over which ordinary linear propagation is forbidden, irrespective of direction. These 

photonic band gaps (PBGs) lend themselves to numerous diversified applications (in linear, nonlinear and 
quantum optics). For instance, PBG structures with line defects can be used for guiding light. Similarly, as 
it has been predicted and confirmed experimentally, photonic crystals allow to modify spontaneous emission 
rate due to the modification of density of quantum states. In particular, it is well known that the density 
of states grows at the edge of the photonic band gap of the PCs. This allows us to predict a higher optical 
gain, but on the other hand a higher level of noise in light generated in PC-lasers operated at a frequency 
near the band gap. 

Photonic crystals are usually viewed as an optical analog of semiconductors that modify the properties of 
light similarly to a microscopic atomic lattice that creates a semiconductor band gap for electrons [3] . It is 
therefore believed that by replacing relatively slow electrons with photons as the carriers of information, the 
speed and bandwidth of advanced communication systems will be dramatically increased, thus revolutionizing 
the telecommunication industry. To employ the high-technology potential of photonic crystals, it is crucially 
important to achieve a dynamical tunability of their band gap [4]. This idea can be realized by changing 
the light intensity in the so-called nonlinear photonic crystals, having a periodic modulation of the nonlinear 
refractive index [5]. Exploration of nonlinear properties of photonic band-gap (PBG) materials could be 
exploited new applications of photonic crystals to devise all-optical signal processing and switching, which 
indicates an effective way to generate tunable band- gap structures by operating entirely with light. 

During the past few years, there has been a great deal of interest in studying propagation of waves inside 
periodic structures. These systems are composites made of inhomogeneous distribution of some material 
periodically embedded in other with different physical properties. Phononic crystals (PCs) [6, 7] are one of 
the examples of these systems. PCs are the extension of the so-called Photonic crystals[8] when elastic and 
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FIG. 1: The motion path of light in one-dimensional function photonic crystals and two-dimensional motion space 

acoustic waves propagate in periodic structures made of materials with different elastic properties. When one 
of these elastic materials is a fluid medium, then PCs are called Sonic Crystals (SC)[9, 10]. For these artificial 
materials, both theoretical and experimental results have shown several interesting physical properties [11]. 
In the homogenization limit [12], it is possible to design acoustic metamaterials that can be used to build 
refractive devices [13]. In the range of wavelengths similar to the periodicity of the PCs, multiple scattering 
process inside the PC leads to the phenomenon of so called Band Gaps (BG), which are required for filtering 
sound [10], trapping sound in defects [14, 15] and for acoustic wave guiding [16]. 

In the present work we present a new kind of function photonic crystals, which refractive index is a function 
of the space position. Unlike conventional PCs, which grow from two materials, A and B, with different 
dielectric constants ea and £b- a periodic layered medium ...A/B/A/B... in case of one-dimensional photonic 
crystals and periodic arrays of cylinders and spheres of the material A embedded in a dielectric matrix B, 
in case of two-dimensional and three-dimensional photonic crystals, respectively. Function PCs may extend 
the concept of PCs, leading likely to some new applications. To exemplify the idea of function PCs, we 
present theoretical calculations of the function photonic crystals band structures and transmissivity. Our 
results indicate that the function photonic crystals behaves more width or more narrow band gap structure 
than the conventional photonic crystals. 

2. The light motion equation in function photonic crystals 

For the function photonic crystals, the crystals refractive index is a periodic function of the space position, 
which can be written as n(z), n(x, z) n(x, y, z) corresponding to one-dimensional, two-dimensional and three- 
dimensional function photonic crystals. In the following, we shall deduce the light motion equations for the 
one-dimensional, two-dimensional and three-dimensional function photonic crystals. Firstly, we give the light 
motion equation in one-dimensional function photonic crystals, and two-dimensional motion space, i.e., the 
refractive index is n — n(z), meanwhile motion path is on xz plane. The incident light wave strikes plane 
interface A point, the curves AB and BC are the path of incident and reflected light respectively, and they 
are shown in FIG. 1. 

The light motion equation can be obtained by Fermat principle, and it is 




(1) 



In the two-dimensional transmission space, the line element ds is 



ds = y/(dx) 2 + {dz) 2 = v 7 ! + z 2 dx, 



(2) 



where z = g§, then Eq. (1) becomes 




(3) 



i.e., 



The Eq. (3) change into 

r° { awz)f+V) Sz + °M*)f+V) Si)dx = 0j (4) 

J A oz 

f B ^l^~pSzdx+ ( B n{z)z(l + z 2 )-^d8z 
J A dz J A 

= [ + z 2 5zdx + n(z)z(l + z 2 )-hz\^ - f d{n{z)z{l + z 2 )-^)5z 

J A dz J A 

= 0. (5) 
The two end points A and B, their variation is zero, i.e., 5z(A) = 5z(B) = 0, and the Eq. (5) is 

f B ,dn(z) r — dn(z) .,,„ i . . zVl + i 2 — z 2 'z(l + i 2 )~5 

/ (— ^Vl + i 2 - — r^z 2 (1 + i 2 )-i -n(z)— ^ '- — )5zdx = {). (6) 

J A dz dz 1 z 

For arbitrary variation 5z, there is 

dn(z) I — dn(z).^. .,. i , , z\J\ + z 2 — z 2 z(l + z 2 )~^ 

simplify Eq. (7), we have 

The Eq. (8) is light motion equation in one-dimensional function photonic crystals and two-dimensional 
motion space. Similarly, we can attain light motion equation in one-dimensional function photonic crystals 
and three-dimensional motion space. It is 

^^y(l + y 2 + z 2 )-n(z){yz-zy) = 0. (9) 
dz 

For the two-dimensional function photonic crystals, light in two-dimensional motion space, the light motion 
equation is 

dz ( ' >i + p dx u ' ( W) 

and the light in three-dimensional motion space, the motion equation is 



dn(x,z) ,.2,-2 



dz 



y(l + y 2 + z 2 ) - n(x, z){yz - zy) = 0. (11) 



For the three-dimensional function photonic crystals, light in two-dimensional motion space, the light motion 
equation is 

Onix^z) On^z). ^ z) _z_ = Q> {u) 

where yo is constant, and the light transmits in three-dimensional motion space, the light motion equation 
is 



(l + y +z )( — z ^ y) - n(x, y, z)(yz - zy) = 0. (13) 
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FIG. 2: The light transmission figure in arbitrary middle medium 
3. The transfer matrix of one-dimensional function photonic crystals 

In this section, we should calculate the transfer matrix of one-dimensional function photonic crystals and 
two-dimensional motion space. In fact, it is the reflection and refraction of light at a plane surface of two 
media with different dielectric properties. The dynamic properties of the electric field and magnetic field are 
contained in the boundary conditions: normal components of D and B are continuous; tangential components 
of E and H are continuous. We consider the electric field perpendicular to the plane of incidence, and the 
coordinate system and symbols as shown in FIG. 2. 

On the two sides of interface I, the tangential components of electric field E and magnetic field H are 
continuous, there arc 



En = Ej = Et: 



E'. 



'0 — Jiil — n>tl -r J^r2 

H = Hi = H tl cos 6l - H' r2 cos 9f, 



(14) 



On the two sides of interface II, the tangential components of electric field E and magnetic field H are 
continuous, and give 



the electric field E t \ is 



En — Ej — E i2 + E r2 

H u = H' I= H a cos 9f - H r2 cos 9f, 



&ti = &tioe v \z=o = ^iioe c w * , 



(15) 



(16) 



and the electric field E i2 is 



E l2 = E tw e i{k '* XB+k '* z) \ z =b = E tlo e l ^ n{bKsin0l ' XB+cos ^ b) . (17) 

Where xa and xb are x component coordinates corresponding to A and B points. 
Now, we calculate the incident angle From Eq. (8), it is straightforward to derive 

dn(z) _ zdz 

-^(z) TTi2' (18) 

then integrate the two sides of Eq. (18) 

r (6) dn(z) _ [ k » zdz 

J n( o) n(z) J k0 1 + z*' [ yj 



to get 



where 



<W - IT**' ,20) 



W 
dz . 



k b =—\ z=b = cot0i, (21) 



dz 

k = —\ z=a = cotel, (22) 



and 



^ = arcsin(-^sin0°), (23) 



where n is air refractive index, and n(0) = n(z)\ z= o- 

By substituting Eqs. (21), (22) and (23) into (20), we attain 

1 + cot 2 0[ _ n 2 (b) 
1 + cot 2 0{ ~ n 2 (0)" 

From Eq. (24), we can find when n(0) = n(b), there is 

o{ = el 



Integrate the two sides of Eq. (18), we can obtain the coordinate component xb 

r n ^ dn{z) _ r k > 
^(o) n{z) J ko 



to get 



since k z > 0, there is 



,n(z) 2 = 1 + fc 2 



i.e., 



\/(l + *o 2 )($}) 2 -l 
Obviously, n(z) > n(0)sin#j. The coordinate ib can be obtained as 



dz 



(24) 



(25) 



/■»<*> rfn(z) = /*■ 

A(0) ™O0 An 1 + ^ 2 ' 



(f^) 2 = fvS. (27) 



^-S^^^O 2 - 1 ' (28) 



dx = dz . (29) 



x B =x A + . (30) 
Jo V(l + ^o 2 )($}) 2 -l 



By substituting Eq. (30) into (17), there is 

uj dz 
E i2 = E tl0 exp[i-n(b)(sm6l \ + smOfxA + cos 0/6)], (31) 

7o V / ( 1 + fc o)(^) 2 -l 

With substituting Eq. (16) into (31), there is 

E i2 = E tl e i5 \ (32) 



where 



rb 

S b = -n(6)(cos 0/ 6 + sin 0f / 

C Jo 



V( 1 + fc o)(^) 2 



and similarly 



(33) 



E 'r2 



E r2 e 



(34) 



Substituting Eqs. (32) and (34) into (14) and (15), and using H = ^Jj^nE, there are 

f Ej = E tl + E r2 e i5 » 



H, = 



<-n(0)E tl cos 6^ - x [^n{Q)E r2 cas6{e iS >> 



and 



En = E tl e iS » + E r2 

H n = .p^n{b)E tl e i5b cos0/ - .[^n{b)E r2 cosflj. 



(35) 



(36) 



From Eq. (35) and (36), we can obtain 
E T 
Hi 



or 



define M matrix 



where 



a„(6) cos fl J ' V " ; 
n(0) cos 



j(Jjf-n(b) cosOf cosSbEu — i sin 5bHii) 



nWcos^ HygHk) co&OlsmSbEn + cos S b H n ), 



Ei 
Hi 



M = 



—i. 



cos <5(, 



f-n(O) cos6» t J sin5 6 



i sin 5h 



n(b) co '' 



n(0) cos ( 



n(fe) cos &l 



T COS (5b 




Hi 



= M 



En 
Hn 



cos Sb 

-i./jf-n(0) cos 0/ sincifc 



fSL n lb) cosS 1 
V w 1 

n(0)cos#f r 

-^T\ COS Of,. 

n(b) cos fft " 



The Eq. (40) is the transfer matrix M of the half period. When n(0) = n(b), there is 



(37) 



(38) 



(39) 



(40) 



(41) 



and the M matrix becomes 



M = 



COS Ob 



sin 5b 



-irjbSmSb cos Sb 



(42) 



where 77^ = y^n(O) cos 0|. In the following, we should study the one-dimensional function Photonic crystals 
at n(0) = n(b). 

4. The structure of one-dimensional function photonic crystals 

In section 3, we attain the M matrix of the half period. We know that the conventional photonic crystals 
is constituted by two different refractive index medium, and the refractive indexes are not continuous on 
the interface of the two medium. We could devise the one-dimensional function photonic crystals structure 
as follows: in the first half period, the medium refractive index is n\{z), and in the second half period, the 
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FIG. 3: The two periods transmission figure of light in function photonic crystals 



medium refractive index is n 2 (z) 7 corresponding thickness b and a, respectively. Their refractive indexes 
satisfy condition n\{b) ^ "2(0), as shown in FIG. 3. We should discuss two kinds of incidence cases: 

(1) The light vertical incidence 

For the vertical incidence, the initial incidence angle 9® = and refraction angle Q\ = 0. From Eq. (41), 
there is 

0! - 0, (43) 

then 

6^=0, el^O, fl t 77/ = 0, ... (44) 

i.e., all incidence angles and refraction angles are zero, and the M matrix is 

M=( C0SS \ ~i si f 6b ). (45) 
\ — iribsmdb cos Ob J 

Where 5b = ^n(b)b and rjb = ^Jj^n(b). In the case, the function photonic crystals becomes conventional 
common photonic crystals, which its refractive indexes n a and rib are constants. 

(2) The light non- vertical incidence 

For light non-vertical incidence, the initial incidence angle 6° 7^ in the first half period I. 
By refraction law, there is 

where n is air refractive indexes, and ni(0) = nx(z)\ z=0 . 

when ni(0) = ri\(b), 0\ = Q\, we obtain the Mb matrix in the first half period I as 

T I cos 5l — -V sin Si \ 

M b = [ ■ I ■ r J [J (47) 

V — if]b sindj coso^ J 



where 



r\{ = rn(O)cos0 t J (48) 
V A*o 



and 

dz j 
- 1 



5{ = -ni(0)[cos6»[ -6 + sin^ 7 / — == 



V( 1 + fc o)(^) 2 
Where fco = cot# 7 . 

In the second half period 77, by refraction law, there is 
when n 2 (0) = n 2 (a), we have 

of 1 = el 1 , 

The M a matrix in the second half period II is obtained 

r ( cos 8" — h sin 6" \ 
M =( vL 1 a \ 

-irfj sin Si 1 cos 5^ /' 



where 



and 



v'J = J-n 2 (0)co S e I t I , 
V Mo 



= -n 2 (0) [cos el 1 ■ a + sin t J/ / 
c Jo 



v / ( i + fc ?)(^) 2 

where fc/ = cot# t . From Eqs. (46) and (50), we have 

sin# 7 = -^-sin0°, 
' n 2 (0) 

and 



cos el 1 



Then the M matrix in the first period is expressed as 



M 1 = Ml ■ M ! a 



In the next, we should calculate the M matrix in the second period. 
By refraction law, there is 



when ni(0) = ni(6), we have 



We can obtain the Mb matrix in the second period 



where 



„ / cosd" 1 — At sin 61" 
Ml 1 = 



V Mo 



si" = %(o) • b + siner f dz ] = si 



when n 2 (0) = n 2 (a), we have 



or = or = o» (es) 



We can obtain the M a matrix in the second half period IV of the second period 

cosS 1 / -^sin<5 
-ir]i v sin Sl v cos^ y 



where 



9 

(62) 



From Eqs. (60)-(62), wc find 

M h u = Mi (63) 

Now, we calculate the M a matrix in the second period. 
By refraction law, there is 



sin Or = ^ sin 0\ u = sin 0° = sin 0\\ (64) 

n 2 (0) n 2 (0) 



cos (5f ^ sin 

M " = ...ivJU/v . (66) 



V Mo 



and 

- -n 2 (0)[cos^ /y a + sin# t /y / cot(ln + arctan(cot 0{ v ))dz\ = S 1 ,, 1 . (68) 

Jo n 2 (0) 



From Eqs. (66)-(68), we find 

Mi 1 = Mi. (69) 

By calculation, we find that all M a matrixes are equal, and all Mj, matrixes are also equal in different period, 
they are 

( cos<5 a — — sin<J a \ 
M a = . .° c , (70) 



and 



where 



-ir/aSinSa cosS a 



M b =[ C0SS \ -i 8 ^ 6 J, (71) 



-irjb sin Sb cos S, 



b 



Va = J^n 2 (0)Jl-^f- S in 2 0l (72) 

V n 2(°) 



'-«i(0)Jl 
Mo V 



n?(0) 



sin 



(73) 



(74) 
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and 

5 b =- ni (0) [cos 6l-b + sin 61 [ dZ =}, (75) 

where fco = cotdj. and fcj = cot Ol 1 . 

Finally, we obtain the M matrix for every period 

M — M b M a — ( 008 \ _ »™* 6 V COSS \ -^ Si ;^V (76) 
\-i?76Sinc)6 cosdf, J \ -ir) a smd a cosd a ) 

For the N-th period, the field vector of up and down En, Hn and En+i , £W+i are satisfied with characteristic 

equation 

En \ — M ( En+1 
Hn J \ Hn+i 



= m n • (77) 



From Eq. (77), we can further obtain the characteristic equation of the N periods photonic crystals, it is 
El ^ M,M 2 ■ ■ ■ M N ( l N+1 ) = M b M a M b M a ■ ■ ■ M b M a ( En+1 



-"(%z)-(oD)(iizy <78 » 

5. The dispersion relation, band gap structure and transmissivity 



With the transfer matrix M (Eq. (76)), we can study the dispersion relation and band gap structure of 
the function photonic crystals. 
From Eqs. (70) and (71), there is 



By Bloch law, we have 



where d = b + a. With Eqs. (79) and (80), there is 



En \ = M b M a ( En +i ) = e~ ikd ( En+1 ) , (81) 



The non-zero solution condition of Eq. (81) is 

det(M b M a - e- lkd ) = 0, (82) 

i.e., 



(cos 5 b cos 5 a — — sin S b sin S a — e lkd )(cos 5 b cos<5 a — — sin S b sin S a — e tkd ) 
r\ b T]a 

+ ( cosJf, sin5 a sin 5 b cos 5 a )(—ir] b s'mS b cos<5 a — irj a cos 5 b s'm5 a ) = 0, (83) 

r\ a Vb 

we resolve the dispersion relation 

coskd = cosS b cos5 a — ~(— + — ) sin S b sin 5 a , (84) 
2 Va Vb 

From Eq. (84), we can study the photonic dispersion relation and band gap structure, and we can obtain 
the transmission coefficient t from Eq. (78) 

t - EtN+1 - ^2 7 (85) 

En Ari + Brj r] N+1 + C + Dr] N+1 ' 

and transmissivity T 

T = t-t* (86) 
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FIG. 4: The picture of three kinds of functions refractive indexes in a period. 



6. Numerical result 

We report in this section our numerical results of transmissivity and dispersion relation. We consider three 
kinds of functions form refractive indexes in a period, 

(1) The first one is sine type function refractive indexes, as 



m(z) = ni(0) + Aisvafz, 0<z<b 
n 2 (z) = n 2 (0) + A 2 sin ^z, 0<z<a 

which are shown in FIG. 4 (a) and (b). 

(2) The second one is upward fold line type function refractive indexes, as 



ni ( ) + 2(m-i)m(0) Z; 0<z< 

m(Q)+ 2(m-l)n l( 0) (b _ z)| |<^< 6j 



and 



2(m 


-l)ni(O) 


2(m 


6 

-l)ni(O) 




6 


2(m- 


-l)ni(O) 


2(m- 


-l)ni(O) 


a 



which are shown in FIG. 4 (c) and (d). 

(3) The third one is downward fold line type function refractive indexes, as 

- ni(0) _ 2n 1 (0)(l-l.lsin^) Z| < Z < § 



and 






6 




2m(0)(l 


— 1.1 sin 






6 




2m(0)(l- 


-1.1 sin t 


D 




a 




2m(0)(l- 


-1.1 sin 6 



ni(Q) _ t2 (6- z ), |<z<6 



z, 0<z<£, 



(87) 



U 2 (0) + ™ Z] 0<z<f, .... 

" 2(Z)= n 2 (0) + ^4^ (a _ z) , § < z < 2 a , ( 8Q ) 



m(*) = < , ..... » . . ' "---2- ( 91 ) 

" ' ; \a-z), %<z<a, 

which are shown in FIG. 4 (e) and (f). 

While ni(0), n 2 (0), m, Ai and A 2 are constants, 6 and a are half period thickness. With Eq. (84), we 
can investigate the dispersion relation and band gap structure, and can resolve transmissivity from Eqs. 
(85) and (86). In FIG. 5, we take sine type function refractive indexes (Eq. (87)), and the parameters are 
09 = f , A 1 = 100, A 2 = 200, m(0) = \/L9, n 2 (0) = y/KE, a = 7A0nm and b = 1260nm. The FIG. 5 (a) 
is the dispersion relation and FIG. 5 (b) is the transmissivity. In the two figures, we can find sine function 
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type photonic crystals has band gap structure. In FIG. 6, we take upward fold line type function refractive 
indexes (Eqs. (88) and (89)), and the parameters are: 9® = ^, ni(0) = -\/h9, n 2 (0) = \/5J5, a = 740nm, 
b = 1260nm and m — 10. The FIG. 6 (a) is the dispersion relation and FIG. 6 (b) is the transmissivity. The 
two figures show the band gap structure when refractive index is upward fold line type function. In FIG. 7, 
there is band gap structure when the refractive index is taken downward fold line type (Eqs. (90) and (91)), 
and the parameters are 9® = J, ni(0) = vT9, n 2 (0) = a = 740nm, b = 1260nm, sm9^ = w "° Q ^ sin 6*°, 

sm9( I — n "° Q ^ sin 9® and no = 1. In FIG. 8, we compare the band gap structures of function photonic 
crystals with the conventional photonic crystals. The FIG. 8 (a) is conventional photonic crystals band 
gap structures, the parameters are: 9® = ^, n(a) = \/F9, n(b) — \/5.5, a = 740nm, and b — 1260nm. 
The FIG. 8 (b) is sine type function photonic crystals band gap structures, the parameters are: 9® = ^, 
A-l = 0.01, A 2 = 0.02, m(0) = \/L9, n 2 (0) = Vhl, a = 740nm and b = 1260nm. The FIG. 8 (c) is also 
sine type function photonic crystals band gap structures, the parameters are: 9® — -|, A\ — 100, A 2 = 200, 
ni(0) = "\/F9, ^2(0) = VKE, a = 740nm and b — 1260nm. The FIG. 8 (d) is upward fold line type function 
photonic crystals band gap structure, and the parameters are: 9® = |, ni(0) = \/L~9, «2(0)) = VbH, 
a = 740nm, b = 1260nm and m = 10. The FIG. 8 (e) is downward fold line type function photonic crystals 
band gap structure, and the parameters are: 9® = |, ni(0) = %/Ll), n 2 (0) = -\Zo\5, a = 740nm, 6 = 1260nm, 
sin 6*^ = k "°q) sinfl°, sin 6^ = n "° ) sin and n = 1. In FIG. 8 (c) and (d), the band gap are more wider 
than the conventional photonic crystals. In FIG. 8 (b) and (e), the band gap are more narrower than the 
conventional photonic crystals. In order to satisfy different application, we can design the different kind of 
function photonic crystals by choosing different refractive index function form. 
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7. Conclusion 

In conclusion, we present a new kind of function photonic crystals, which refractive index is a function of 
space position. Unlike conventional PCs, which structure grow from two materials A and B, with different 
dielectric constants sa and sb- Based on Fermat principle, we achieve the motion equations of light in 
one-dimensional, two-dimensional and three-dimensional function photonic crystals. For one-dimensional 
function photonic crystals, we investigate the dispersion relation, band gap structure and transmissivity, 
and compare them with conventional photonic crystals. By choosing different refractive index distribution 
function n(z), we can obtain more wider or more narrower band gap structure than conventional photonic 
crystals. Due to the function photonic crystals has more wider or more narrower band gap structure 
than conventional photonic crystals, we think the function photonic crystals should has more extensive 
application foreground. 
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FIG. 6: The dispersion relation, band gap structure and transmissivity for upward fold line function refractive indexes 
(Eqs. (88)and (89)). 
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FIG. 8: Compare the band gap structures of different refractive index distribution function n(z) function photonic 
crystals with conventional photonic crystals. 



